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The consumption of insects by apes has previously been reported based on direct observations and/or trail 
signs in feces. However, DNA-based diet analyses may have the potential to reveal trophic links for these 
wild species. Herein, we analyzed the insect-diet diversity of 9 feces obtained from three species of African 
great apes, gorilla (Gorilla gorilla gorilla), chimpanzee (Pan troglodytes) and bonobo (Pan paniscus), using 
two mitochondrial amplifications for arthropods. A total of 1056 clones were sequenced for Cyt-b and COI 
gene libraries, which contained 50 and 56 operational taxonomic units (OTUs), respectively. BLAST 
research revealed that the OTUs belonged to 32 families from 5 orders (Diptera, Isoptera, Lepidoptera, 
Coleoptera, and Orthoptera) . While ants were not detected by this method, the consumption of flies, beetles, 
moths, mosquitoes and termites was evident in these samples. Our findings indicate that molecular 
techniques can be used to analyze insect food items in wild animals. 

I nvestigating the diets of primate populations elucidates their behavioral ecology and evolution and clarifies 
I dietary differences among the same species in different habitats and among different species in shared 
I habitats 1 . The consumption of insects is widespread among nonhuman primates, and their entomophagy 
depends on their body size. Small primates such as galagos (Galago crassicaudatus and Galago senegalensis), 
pottos (Perodicticus potto), and tarsiers (Tarsius spectrum) are obligate insect eaters; their diet is composed mainly 
of insects belonging to the orders Lepidoptera, Orthoptera and Hymenoptera 24 . The diets of medium-bodied 
primates such as red colobus monkeys (Procolobus tephrosceles) and blue monkey (Cercopithecus mitis) consist 
mainly of young leaves, flowers, and unripe fruit; they are also known to eat insects less frequently than small 
primates 5,6 . Although the large-bodied great apes have a high frugivory 711 and folivory component to their 
diet 12,13 , there is also an insectivorous component in the diet of gorilla, chimpanzee, bonobo, orangutan and 
gibbon populations 14 20 . 

Insects comprise 0.1-1% of the fresh weight of the daily food ingested by both gorillas and chimpanzees 21 , but 
they may provide disproportionate nutritional benefits due to their role as protein sources. Insects are particularly 
valuable because they provide certain amino acids, vitamins (such as B12) and minerals (including iron (Fe) and 
manganese) that may be absent in plant foods 17,22 24 . 

The insect prey items vary among gorilla species living in different habitats. In mountain gorillas (Gorilla 
beringei beringei), insectivory is uncommon 18,25,26 , but they feed on ants (Dorillus spp.) to a variable degree 18,26 . In 
contrast, Grauer's gorillas (Gorilla beringei graueri) consume ants more regularly than mountain gorillas, perhaps 
because the insects are more prevalent 27 . Furthermore, western gorillas (Gorilla gorilla gorilla) feed on insects, 
indicating that insects are an important aspect of the diet of these gorillas 10,14,16,17,28,29 . 

There is considerable temporal variation in insect consumption among different chimpanzee populations. 
Both forest-inhabiting chimpanzees in central Africa and savannah-dwelling chimpanzees in western Africa 
consume insects, such as termites, throughout the year 21,30 32 . In East Africa and at Mt. Assirik in Senegal, 
chimpanzees feed on termites only seasonally 33 . There is not much variation in insect consumption by different 
chimpanzee groups. Generally, chimpanzees eat ants, termites, bees, wasps, caterpillars and beetle grubs 29,34,35 . In 
contrast, arthropods comprise a small fraction of the diet of wild bonobos. The most important non-plant food 
sources for bonobos are invertebrates, including larvae, termites, ants, bees, earthworms and millipedes 36-39 . 
Insectivory is also reported for orangutans; especially during times of fruit scarcity 40 , orangutans start to prey 
on insects including termites, ants, bees, gall wasps, crickets, caterpillars and bush crickets 19,41,42 .Wild Sumatran 
orangutans (Pongo pygmaeus abelii) make and use tools (twigs, sticks and branches) for extractive insect 
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foraging 41 . Very little data is available on the diet of gibbons, although 
fruit is their main food source, which is supplemented with some 
arthropods including caterpillars and termites 20,43 . 

The studies cited above are based either on direct observations of 
feeding apes or morphological classification of insect remains in ape 
feces. Another avenue for investigation of insectivory is the study of 
feeding remains left by apes (e.g. to identify insect species eaten from 
a recently disturbed ant nest or encountered termite mound). 
Although these classical methods can provide direct records for both 
qualitative and quantitative behavior (i.e. insect prey choice and 
frequency of insect consumption) and they are inexpensive. 
However, direct observations of insects consumed by apes are limited 
in each individual trial by difficulties in tracking apes in their native 
habitat; this approach can be time-consuming and therefore 
impractical for studying insectivory in non-human primates 44 . 
Furthermore, the damage to insect bodies and distortion of their 
morphology by the mechanical action of the digestive tract creates 
major challenges for direct analyses of insects in animal feces 44 . 
Alternative molecular methods may provide better identification 
and global detection of insects in primates feces compared with these 
classical methods and may also elucidate the feeding behavior of 
primates 45,46 . 

The first PCR-based method for detecting animal prey in primate 
feces was used by Hofreiter et al. 46 , who detected fragments of verte- 
brate mitochondrial DNA in fecal samples of gorillas and bonobos. 
In a metagenomic study with primers targeting the arthropod cyto- 
chrome b (Cyt-b) gene, Pickett et al. 45 investigated the insect diets of 
six sympatric New World monkeys living in a western Amazon rain- 
forest. To date, molecular analyses of great ape insectivores using 
high throughput sequencing have not been reported. However, using 
this technique to target the cytochrome c oxidase I (COI) fragment is 
effective for studying the arthropod diet of bats and, thus, could be 
applicable to other insectivores 47,48 . 

In this study, we evaluated the insect-diet diversity of three African 
great apes gorilla, chimpanzee and bonobo by analyzing fecal sam- 
ples (N = 9) with DNA-barcoding primers targeting the Cyt-b 49 and 
COP° genes in arthropod mitochondrial genomes. 

Results 

A total of 1056 clones were collected from the cloning libraries gen- 
erated by using both primer sets. First, 24 clones per sample and per 
amplification were sequenced and analyzed. No additional opera- 
tional taxonomic units (OTUs) 51 were found when we sequenced 
24 additional clones from each sample. Moreover, in one sample 
(Gl), sequencing 144 clones did not increase the number of OTUs 
that were obtained from the first 24 sequenced clones. Finally, after 
examining the 1056 clones, 1006 clones presented good quality 
sequences that did not contain termination codons within the trans- 
lated sequences. Analyses of these 1006 clones resulted in the detec- 
tion of 106 different insect OTUs in the 9 fecal samples from bonobos 
(Bl, B2 and B3), chimpanzees (CI, C2 and C3) and gorillas (Gl, G2 
and G3) with both Cyt-b and COI genes (Fig. 1 and 2). 

Analysis of the Cyt-b gene revealed 18 OTUs from insects belong- 
ing to 12 families in 4 orders (Isoptera, Diptera, Lepidoptera and 
Coleoptera) in the fecal samples from bonobos (Fig. 1 and 
Table 1). Nineteen OTUs belonging to 10 families from Diptera 
and Coleoptera and 16 OTUs belonging to 1 1 families from 5 orders 
(Diptera, Isoptera, Lepidoptera, Coleoptera and Orthoptera) were 
identified in the feces of chimpanzees and gorillas, respectively, with 
the same primer set (Fig. 1 and Table 1). 

The COI bonobo clone libraries yielded 23 OTUs belonging to 12 
families within the orders Diptera and Lepidoptera (Fig. 2 and 
Table 2). Similarly, 19 and 21 OTUs from 12 and 10 families within 
the orders Diptera, Lepidoptera and Coleoptera were retrieved from 
the COI chimpanzee and gorilla clone libraries, respectively (Fig. 2 
and Table 2). 



Taking the PCR results together, we detected a variety of arth- 
ropod OTUs in the fecal samples from these three African great apes, 
including fruit flies, moths, beetles, butterflies, mosquitoes and ter- 
mites (Fig. 1 and 2 and Tables 1 and 2). A total of 32 families from 5 
orders were present in at least one of the African great apes studied 
(Fig. 3). Eight of these families (Carabidae, Chrysomelidae and 
Staphylinidae from the order Coleoptera and Cecidomyiidae, 
Muscidae, Psychodidae, Sciaridae and Sepsidae from the order 
Diptera) were commonly detected in feces from the 3 African great 
apes species (Fig. 3). Three types of OTUs belonging to the order 
Isoptera (termites) were only identified in gorilla and bonobo feces 
(Gl, G2 and B2), as shown in Fig. 1 and 3. However, no OTUs 
assigned to the families Formicidae, Apidae and Tabanidae were 
identified in our samples (Fig. 3). 

Discussion 

Investigating the diets of animals by applying molecular methods to 
fecal samples is useful for the study of wild animals that are difficult 
to observe 52 54 . Very few studies have examined the insect-diet 
diversity eaten by primates through molecular approaches. 
Hofreiter et al. 46 investigated the presence of vertebrate prey DNA 
in bonobo and gorilla feces through PCR-based methods, and more 
recently, Pickett et al. 45 evaluated insect diversity in fecal samples 
from new world monkeys using a single arthropod primer set. 
Thus, our study is the first to analyze insect diversity in fecal samples 
from African great apes using two primers (targeting the Cyt-b and 
COI genes in the mitochondrial genome of arthropods). These pri- 
mers have been used to successfully analyze the arthropod prey of 
tiger beetles and bats 49,50 . In this study, we detected 106 insect OTUs 
(41, 38 and 37 for bonobos, chimpanzees and gorillas, respectively) 
using these primers (Fig. 1 and 2, Tables 1 and 2). Compared with 
behavioral observations and/or analysis of trail signs in ape feces, we 
found many previously unknown insect families that are consumed 
by African great apes (Fig. 3). Many insects, such as species in the 
orders Coleoptera and Lepidoptera or caterpillars detected in this 
study (Tables 1 and 2) have strong associations with plants 55,56 . Some 
of these insect species, such as member of family Chrysomelidae, feed 
on different plant parts 57 . Consequently, they could be eaten incid- 
entally (secondary predation) when African great apes feed on 
plants. Thus, one advantage of using a molecular approach to exam- 
ine insects consumed is the inclusion of those consumed via second- 
ary predation. These species may not have been otherwise detected 
through classical approaches, but they are still components of the diet 
that have nutritional value. These indirectly eaten insects may also 
contribute to an understanding of the feeding ecology and foraging 
strategy of a species. Interestingly, different OTUs from termites 
were identified in both gorilla and bonobo feces (Fig. 1 and 
Table 1). Termites are commonly preyed on by African great apes 21,33 
(Fig. 3). Macrotermes species, which are the largest termites in Africa 
and live in different habitat types, are an important insect meal for 
both western gorillas and chimpanzees residing in southeastern 
Cameroon 21 (Fig. 3). Despite the wide global distribution of some 
of the insect species found in our work, such as those belonging to 
Drosophilidae, no data on the occurrence of the remaining species in 
Cameroon and DRC are available. 

We were unable to amplify any of Formicidae species in this study 
(Fig. 3 and Tables 1 and 2) using both arthropod primers, and the 
macroscopic examination of fecal samples did not detect any visible 
hard parts (e.g. chitinous exoskeleton or jointed appendages) of these 
insects. Their absence could be due to degradation of insect DNA 
during its journey through the digestive tube, the small sample size or 
insufficient mitochondrial DNA caused by bias in the DNA extrac- 
tion, amplification and cloning process. 

While DNA-based fecal analyses are useful for evaluating the diets 
of wild animals, several challenges and limitations should be noted 
when employing this method. Although, no hard insect parts were 
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Figure 1 | Insect OTUs detected in African great ape fecal samples using CB3/CB4 primers targeting the Cyt-b gene. The neighbor-joining tree was 
constructed using MEGA 5 and supported by 500 bootstrap replicates. The insect species listed with accession numbers are the closest matching 
sequences from GenBank. Gl, G2 and G3: fecal samples from gorillas; CI, C2 and C3: chimpanzees; Bl, B2 and B3: bonobos. 
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Figure 2 | Insect OTUs detected in African great fecal samples using ZBJ-ArtFlc/ZBJ-ArtR2c primers targeting the COI gene. The neighbor- 
joining tree was constructed using MEGA 5 and supported by 500 bootstrap replicates. The insect species listed with accession numbers are the closest 
matching sequences from GenBank. Gl, G2 and G3: the fecal samples from gorillas; CI, C2 and C3: chimpanzees; Bl, B2 and B3: bonobos. 
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1 r\ 1 rA nn i-> —r i-\ i-\ 
V . 1 1 1 r T 1 UU 1 1 zee 


C3 


17 f 1 Al 
1 / \ 1 OJ 


1 nlfinntflrn 
V^UIcUUIcl U 


ll\/ti c/~ 1 n n 0 
Lsy 1 1 1 UUC 


Root oc (\A/ntor noot ocl 
Dccllco IVVUIel Uccllcol 


Chimpanzee 




90 (91 


Coleoptera 


Scarabaeidae 


Beetles (Scarab beetles) 


1 1 1 1 1 UU 1 iz.ee 


C3 


/in (91 

4U [ZJ 


C nlfinntflrn 
^-UlcUUIcl U 


llwti c/~ 1 n n 0 
L'y 1 1 5CI UUe 


Root oc (\A/ntor noot ocl 
Dccllco IVVUIel Uccllcol 


x_. 1 1 1 1 1 1 uu 1 iz.ee 


C3 


9fl Ml 
ZU 


1) 1 ntd rn 
\s 1 uiei u 


t n 1 rn n n m 1 n n 0 
x_. 1 1 1 1 u 1 iui 1 1 1 uue 


F 10c Innnnitinn minnocl 
1 neo 1 1 iui iui 1 1 1 iu iinuycoi 


in 1 m n 700 
x_. 1 1 1 1 1 1 uu 1 iz.ee 


C3 


zll (91 


1) 1 nto rn 

\-S 1 Ul CI u 


10/" innmvi inno 
v cl luuiiiy iiuuc 


F 10c li^nll nnntel 
1 neo |Vw/un yiiuioi 


Gorilla 


\j? 1 


1 1 99\ 


1 so pte ra 


Termitidae 


Termite 


Gorilla 


\j 1 


9 1 (91 


1 so pte ra 


Termitidae 


Termite 


Gorilla 


\j 1 


A /Ol 


D i pte ra 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 




AO (91 
4Z 


Coleoptera 


Carabidae 


Beetles 




G2 


y| 9 (91 


(n onntorn 

V UICUUI Cl u 


l nr\/cnmol inno 
x_. 1 11 y oui 1 tei iuuc 


Root oc M onr noot oc 1 
ueeneo ilcui ueeneoi 


Gorilla 


G2 


1 O f Al 


(Irmnntorn 

\ ' 1 II IUUI Cl u 


Acrininno 

AAl_ 1 IUIUUC 


Nrncc nnnnorc 

\J\ Uiol IUUUCI 0 


Gorilla 


G2 


21 (4) 


1) i ntp rn 
1 uiei u 


Son? irlno 


F ioc iFncinn r iocl 

1 1 1 CO ILI IOIUI 1 1 1 1 CO J 


„ ... 
Oonlla 




9A (zll 


Di ptera 


Muscidae 


Flies (House flies) 


Gorilla 


G2 


46 (2) 


Di ptera 


Muscidae 


Flies (House flies) 


Gorilla 


G2 


9/1 Ml 
Z4 14) 


Lepidoptera 


Arctiidae 


Moths (tiger moths) 


Gorilla 


G2 


5 (8) 


1 so ptera 


Rhinotermitidae 


Termite 


Gorilla 


G3 


6 (9) 


Di ptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


G3 


13 (9) 


Di ptera 


Cecidomyiidae 


Flies (Gall gnats) 


Gorilla 


G3 


14(6) 


Di ptera 


Sepsidae 


Flies (Ensign flies) 


Gorilla 


G3 


49 (3) 


Di ptera 


Fanniidae 


Flies (True flies) 


Gorilla 


G3 


9 (21) 


Coleoptera 


Carabidae 


Beetles 



Gl , G2 and G3: fecal samples from gorillas; Cl , C2 and C3: chimpanzees; Bl , B2 and B3: bonobos. 



visible in feces prior to analysis, fecal sample contamination is pos- 
sible and, unfortunately, difficult to investigate especially for insects 
that burrow into feces. However, we minimized this limitation by 
using the inner part of the fecal mass for DNA extraction, which 
reduced contamination with soil organisms, flies and their eggs 
deposition. 

DNA extraction, PCR amplification and cloning biases were also 
limitations in this study. The amount of feces used for DNA extrac- 
tions represented only a small fraction of the total feces bulk. Low 
DNA concentration, species with low frequency in samples and the 
limited number of tested clones also may have prevented the detec- 
tion of certain prey items. Moreover, the primers used may not have 



amplified the mitochondrial sequences of some insect species, so 
these sequences were lost. Another challenge with molecular diet 
analyses using DNA from feces is that one cannot determine whether 
the insects present in the guts of African great apes were the result of 
primary or secondary predation or contamination from the edible 
parts of plants. Although the primer sets we used amplified a wide 
range of arthropods, the incompleteness of insect sequence data in 
GenBank prevented us from assigning the recovered OTUs at genus 
or species level. For this reason, the taxonomic identifications were 
done at family level according to the E-values of BLAST results. 
Exhaustive studies using molecular techniques along with con- 
ventional morphological taxonomy are needed to enrich the public 
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Table 2 | Insect OTUs detected in 9 fecal samples from wild African great apes using CO/ targeting primers 



Ape 


Sample 


OTU (No 


of dones) 


Order 


Family 


Common name 


Bonobo 


Bl 


1 


(31) 


Diptera 


Tephritidae 


Flies (Fruit flies) 


Bonobo 


Bl 


2 


(4) 


Diptera 


Psychodidae 


Flies (drain flies) 


Bonobo 


Bl 


3 


(7) 


Diptera 


Tephritidae 


Flies (Fruit flies) 


Bonobo 


Bl 


4 


(2) 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Bonobo 


Bl 


5 


(2) 


Diptera 


Tephritidae 


Flies (Fruit flies) 


Bonobo 


Bl 


6 


(2) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Bonobo 


B2 


7 


(38] 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Bonobo 


B2 


8 


(2) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Bonobo 


B2 


9 


(4) 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Bonobo 


B2 


1 1 


(2) 


Diptera 


Chironomidae 


Flies (nonbiting midges) 


Bonobo 


B2 


10 


(2) 


Lepidoptera 


Noctuidae 


LA 1 1 / /— \ 1 , , 1 \ 

Moths (Owlet moths) 


Bonobo 


B3 


1 2 


(16) 


Diptera 


Muscidae 


Flies (House flies) 


Bonobo 


B3 


14 


(2) 


Diptera 


Culicidae 


Mosquitoes 


Bonobo 


B3 


1 6 


(2) 


Diptera 


Curtonotidae 


Flies (Quasimodo flies) 


Bonobo 


B3 


1 8 


(2) 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Bonobo 


B3 


1 9 


(2) 


Diptera 


Muscidae 


Flies (House flies) 


Bonobo 


B3 


22 


(2) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Bonobo 


B3 


1 3 


(4) 


Lepidoptera 


Sphingidae 


k i i 1 /II 1 ill 

Moths (Hawk moths) 


Bonobo 


B3 


15 


(6) 


Lepidoptera 


Nymphalidae 


Butterflies (Brush-footed 










butterflies) 


Bonobo 


B3 


1 7 


(4) 


Lepidoptera 


Noctuidae 


i 1 ,| 1/ — > | , ■ 1 I 

Moths (Owlet moths) 


Bonobo 


B3 


20 


(4) 


Lepidoptera 


Pyralidae 


Moths (snout moths) 


Bonobo 


B3 


21 


(2) 


Lepidoptera 


Noctuidae 


Moths (Owlet moths) 


Bonobo 


B3 


23 


(2) 


Lepidoptera 


Noctuidae 


Moths (Owlet moths) 


Chimpanzee 


CI 


1 2 


(2) 


Diptera 


Muscidae 


Flies (House flies) 


Chimpanzee 


CI 


16 


(4) 


Diptera 


Curtonotidae 


Flies (Quasimodo flies) 


Chimpanzee 


CI 


24 


(8) 


Diptera 


Psychodidae 


Flies (drain flies) 


Chimpanzee 


CI 


25 


(8) 


Diptera 


Sciaridae 


Hies (Dark-winged tungus gnats) 


Chimpanzee 


CI 


26 


(4) 


Diptera 


Curtonotidae 


Flies (Quasimodo flies) 


Chimpanzee 


CI 


28 


(2) 


Diptera 


Lecidomyiidae 


Flies (Gall gnats) 


Chimpanzee 


CI 


27 


(20) 


Lepidoptera 


Geometridae 


Moths (geometer moths) 


Chimpanzee 


C2 


29 


(26) 


Diptera 


Tachinidae 


Flies (True flies) 


Chimpanzee 


C2 


30 


(6) 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Chimpanzee 


C2 


32 


(2) 


Diptera 


Tachinidae 


Flies (True flies) 


Chimpanzee 


C2 


35 


(2) 


Diptera 


Tachinidae 


Flies (True flies) 


Chimpanzee 


C2 


3 1 


(2) 


Lepidoptera 


Geometridae 


Moths (geometer moths) 


Chimpanzee 


C2 


33 


(4) 


Lepidoptera 


Cosmopterigidae 


Moths (Cosmet moths) 


Chimpanzee 


C2 


34 


(4) 


Lepidoptera 


Erebidae 


Moths 


Chimpanzee 


C2 


36 


(2) 


Lepidoptera 


Cosmopterigidae 


Moths (Cosmet moths) 


Chimpanzee 


C3 


1 0 


(5) 


Lepidoptera 


Noctuidae 


i i .1 / V 1 . ill 

Moths (Owlet moths) 


Chimpanzee 


C3 


16 


(28) 


Diptera 


Curtonotidae 


Flies (Quasimodo flies) 


Chimpanzee 


C3 


37 


(11) 


Diptera 


Simuliidae 


r— i ■ mi l ri ■ \ 

Flies (Black flies) 


Chimpanzee 


C3 


38 


(4) 


Coleoptera 


Carabidae 


Beetles 


Gorilla 


Gl 


6 


(1) 


Diptera 


Muscidae 


Flies (House flies) 


Gorilla 


Gl 


19 


(4) 


Diptera 


Muscidae 


Flies (House flies) 


Gorilla 


Gl 


39 


(74) 


Diptera 


Muscidae 


Flies (House flies) 


Gorilla 


Gl 


40 


(57) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


Gl 


41 


(3) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


Gl 


42 


(3) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


Gl 


43 


(2) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


G2 


44 


(13) 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Gorilla 


G2 


46 


(6) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


G2 


47 


(2) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


G2 


48 


(10) 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Gorilla 


G2 


49 


(6) 


Coleoptera 


Staphylinidae 


n 1 1 /n 1 1 1 \ 

Beetles (Kove beetles) 


Gorilla 


G2 


45 


(11) 


Lepidoptera 


Noctuidae 


Moths (Owlet moths) 




G3 


44 




1) i r~itd rn 


V- CL lUUMIyllVJVJC 


F hoc i(-tn fi n fi t c 1 
1 Ilea lx_7UII yjllUlol 


Gorilla 


G3 


50 


(6) 


Diptera 


Sepsidae 


Flies (Ensign flies) 


Gorilla 


G3 


51 


(4) 


Diptera 


Drosophilidae 


Flies (Fruit flies) 


Gorilla 


G3 


52 


(4) 


Diptera 


Syrphidae 


Flies (Hoverflies) 


Gorilla 


G3 


54 


(4) 


Diptera 


Psychodidae 


Flies (drain flies) 


Gorilla 


G3 


55 


(2), 


Diptera 


Cecidomyiidae 


Flies (Gall gnats) 


Gorilla 


G3 


56 


(2) 


Diptera 


Tephritidae 


Flies (Fruit flies) 


Gorilla 


G3 


53 


(2) 


Lepidoptera 


Riodinidae 


Butterflies (metalmarks) 


G 1 , G2 and G3 : feca 


samples from 


gorillas; CI , C2 ar 


d C3: chimpanzees; Bl , B2 and B3: bonobos. 
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Figure 3 | Summary of insect consumption by three species of African 
great apes. Orders and families of known insect prey identified through 
DNA analysis, direct observation or morphological analysis of prey 
remains in feces. Arthropod families obtained from Cyt-b and COI 
sequences and from different individuals of the same species were pooled 
together. Green cells indicate insect families that were identified as diet 
items in this study. Red cells indicate insect families that were not detected 
in this study. Blue cells indicate insect families that were detected in this 
study by DNA and in other studies through classical methods. 

databases with complete arthropod mitochondrial genome sequ- 
ences, which will identify insects at the species level. Finally, prey 
sequence counts recovered using traditional cloning/sequencing or 
pyrosequencing 58 cannot quantify insect consumption or provide 
insights to insect choice by African great apes. Thus, combining 
molecular tools with classical observation methods provides a more 
complete picture of the complex diet of wild animals, including 
primates. 

In conclusion, molecular fecal analysis is an emerging method that 
can provide valuable insight into the insectivory aspect of primate 
diet and can expand the list of insect foods ingested by a study 
population. The high number of insect OTUs detected in fecal 
samples using this strategy confirms the importance of molecular 



applications for analyzing direct and indirect insect consumption in 
African great apes. For detailed descriptions of the occurrence and 
abundance of insects in African great ape feces using next generation 
sequencing, additional studies using larger numbers of samples col- 
lected in different seasons are needed. This approach will clarify 
intra- and inter-species dietary variation in wild primates. 

Methods 

Source of the fecal samples. A total of 9 fecal samples were used in this study. Three 
fecal samples from wild gorillas (G. g. gorilla) were collected at a forest site near 
Messok (sample Gl) and at a forest site near Mambele (samples G2, G3) (Messokand 
Mambele are two towns located southeast and east of Cameroon, respectively). Three 
fecal samples from chimpanzees (P. troglodytes) were collected at a forest site near 
Mambele (samples Cl, C2) and at a forest site near Messok (sample C3), and three 
fecal samples (Bl, B2 and B3) from bonobos (P. paniscus) were collected at the 
Lomako-Yokokala faunal reserve in the Democratic Republic of Congo (DRC). The 
sample collection protocol was described previously 59 . Briefly, samples were collected 
near night nests or feeding sites. For all samples, the time, date and GPS position were 
recorded, and the species of origin was determined according to prints, collection sites 
and morphological aspects of the samples as well as by amplifying a 386 bp mtDNA 
fragment spanning the 12SrRNA gene (using primers 12S-L1091 5'-AAAAAGCT- 
TCAAACTGGGATTAGATACCCCACTAT-3 ' and 12S-H1478 5'-TGACTGC- 
AGAGGGTGACGGGCGGTGTGT-3'). The time between defection and collection 
was estimated at <24 h, according to the physical texture of the samples. Collected 
samples were saved at base camps at ambient temperature in KNAlater (Ambion, 
Austin, TX) for less than 3 weeks and transported to a central laboratory for storage at 
— 80°C. No experimentation was conducted on these animals. The collection of fecal 
samples from the soil was approved by the Ministry of Scientific Research and 
Innovation of Cameroon and DRC. No other permit was required, as this research 
was non-invasive work, and the collection of the samples did not disrupt wild fauna. 

DNA extraction. The outer layer of fecal bulk was peeled carefully with a sterile 
scalpel and polystyrene tweezers (approximately 1-2 mm was removed). The inner 
part of the fecal bulk was used for extraction to avoid a possible contamination with 
soil organisms and/or the risk of egg deposition by some flies, which would result in 
the amplification of non-prey organisms. DNA was extracted using a modified 
version of the Qiagen fecal procedure (QIAamp DNA Tissue Kit, Qiagen Inc., 
Germany) 60 . A 200-mg aliquot of each fecal sample was placed in a 2-ml tube 
containing 200 mg of a mixture of 0.1-, 0.5-, and 2-mm zirconium beads and 1.5 ml 
of ASL buffer (Qiagen). The sample was bead-beaten at 3200 rpm for 90 seconds, 
followed by heating at 95°C for 10 minutes. The final pellet was suspended in 180 ul 
of tissue lysis buffer and incubated with proteinase K for 2 hours at 55°C. The 
manufacturer's recommendations were followed for the purification and elution of 
DNA. 

Primer selection and genomic amplification. To minimize primer bias, two 
different primers targeting Cyt-b and COI genes in the mitochondrial genome of 
arthropods were used 49,50 . The first PCR was previously used to study the prey of tiger 
beetles after whole DNA extraction from predatory specimens, which provided lower 
identity sequences compared with the GenBank database but sufficiently high scores 
to identify the order, family and occasionally the genus 49 . The second PCR was used to 
analyze arthropod prey in bat feces, enabling the authors to identify 37 prey taxa from 
15 fecal samples (identification at the species level in 72% of the analyzed clones) 50 . 
The primer sets ZBJ-ArtFlc ( AG AT ATTGG AAC WTT AT ATTTT ATTTTTGG ) / 
ZBJ-ArtR2c ( WACT A ATC A ATT WCC AA ATCCTCC ) 49 and CB3 (GAGGAGCAA- 
CTGTAATTACTAA)/CB4 (AAAAGAAARTATCATTCAGGTTGAAT) 50 
amplified 157 and 410 bp segments of mitochondrial DNA, respectively. The PCR 
reaction mixture (final volume 50 ul) contained 5 ul of dNTPs (2 mM of each 
nucleotide), 5 ul of 10X DNA polymerase buffer (QIAGEN, Courtaboeuf, France), 
1 ul of MgCl 2 (25 mM), 0.25 ul of HotStarTaq DNA polymerase (5 U) (QIAGEN, 
Courtaboeuf, France), 1 ul of each primer (10 pmol/ul), and 5 ul of extracted DNA. 
The PCR was performed with an initial denaturation at 95°C for 15 minutes, followed 
by 40 cycles of 95"C for 45 seconds, an annealing temperature specific for the primers 
used (57°C for ZBJ-ArtFlc/ZBJ-ArtR2c and 46°C for CB3/CB4) for 30 seconds, 72°C 
for 1 minute, and a final extension at 72°C for 5 minutes. The PCR products were 
analyzed using agarose gel electrophoresis (1.5%) and visualized with ethidium 
bromide staining. Positive PCR products were subsequently purified using the 
NucleoFast® 96 PCR Kit (MACHEREY-NAGEL, Hoerdt, France) according to the 
manufacturer's instructions. 

Cloning procedures and insert amplification. The purified PCR products were 
cloned with the pGEM®-T Easy Vector System 2 Kit (Promega, Madison, USA), as 
recommended by the manufacturer. All white colonies were analyzed by PCR with 
the M13d (5 ' -GTAAAACGACGGCCAG) and M13r (5 ' -CAGGAAACAGCTA- 
TGAC) primers, as described previously 60 . 

Sequencing and data analysis. Correct sizes of purified PCR-M13 inserts were 
sequenced in both directions using the Big Dye® Terminator vl.l Cycle Sequencing 
Kit (Applied Biosystems, Courtaboeuf, France). The M13d and M13r primers were 
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used for sequencing. The sequencing products were run on an ABI PRISM 3130 
automated sequencer (Applied Biosystems, Foster city, CA, USA). 

The sequences were corrected with the CodonCode Aligner software (4.0.1). 
Sequence similarity was determined through multiple alignment software using 
ClustalW2. Sequences with base-pair mismatches, insertions or deletions leading to 
stop codons in every possible reading frame were excluded from further analysis. A 
cloned sequence was designated an "operational taxonomic unit" (OTU) 51 if it had 
more than 2% sequence divergence from all other recovered sequences across 
approximately 400 bp and 157 bp of Cyt-B and COI, respectively 45 . Phylogenetic 
analyses (neighbor-joining tree) were performed using MEGA 5, supported by 500 
bootstrap replicates. 

To make taxonomy assignments, the OTUs were compared with a BLAST database 
of pre-assigned sequences in GenBank (available at the National Center for 
Biotechnology Information website: http://www.ncbi.nlm.nih.gov/), and the closest 
available matching hit with the lowest E-value was assigned as the taxonomic iden- 
tification for each OTU. As many arthropod species and genera shared or had very 
close E-values, OTUs were assigned to the family level. 

Nucleotide sequence accession numbers. All sequences generated from the primer 
set CB3/CB4 were deposited in the GenBank database with the accession numbers 
JX500757 to JX501234. The sequences obtained from the primer set ZBJ-ArtFlc/ZBJ- 
ArtR2c were not deposited because sequences shorter than 200 bp are not accepted in 
the GenBank database. 
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